Measurements of the continuum -R u ds and R values in e + e annihilation in the energy 

region between 3.650 and 3.872 GeV 
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We report measurents of the continuum R u d B near the center-of-mass energy of 3.70 GeV, the 
-Ruds(c)+»/>(3770) (s) and the -Rhad(s) values in e + e~ annihilation at 68 energy points in the energy 
region between 3.650 and 3.872 GeV with the BES-II detector at the BEPC Collodier. We obtain 
the f? u ds for the continuum light hadron (containing u, d and s quarks) production near the DD 
threshold to be R udB = 2.141 ± 0.025 ± 0.085. 
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Precision measurement of the cross section c° ad (s) for 
single photon exchange process in e + e~ annihilation at 
the cm. (center-of-mass) energy y/s, which does not in- 
clude the initial state radiative and vacuum polarization 



turbative QCD (pQCD) calculation on the cross section 
and in calculation of the photon vacuum polarization. 
This cross section is defined as the lowest order cross 
section for inclusive hadron production in the paper. It 



2 



is often represented by a R(s) ratio 



R(s) 



'had 



(s) 



(i) 



where <r° + (s) = 47ra 2 (0)/3s is the lowest order QED 
cross section for muon pair production. The cr{J ad (s) still 
includes the final state photon radiation and all QCD cor- 
rections. It can be extracted from the observed cross sec- 
tion crjjj(s) for inclusive hadron production accounting 
for the initial state radiative corrections and the vacuum 
polarization corrections (see below). 

In the lower energy region, precise measurements of 
the iiuds(s), which is defined as the R(s) ratio for the 
continuum light hadron (containing u, d and s quarks) 
production in e + e~ annihilation in this paper, can be 
used to test the validity of the pQCD calculation [l],0 in 
this energy region. Moreover the i?had(s) values includ- 
ing the contributions from both the continuum hadrons 
and all 1 resonances at all energies are needed to calcu- 
late the effects of vacuum polarization on the parameters 
of the Standard Model, such as the quantities a(M§), the 
QED running coupling constant evaluated at the mass of 
Z , and = (g — 2)/2, the anomalous magnetic mo- 
ment of the muon [H, 0, [H . A large part of uncertainty in 
the calculation arises from the uncertainties in the mea- 
sured inclusive hadronic cross sections in the open charm 
threshold region, in which many resonances overlap. To 
get credible measurements of the i? u ds(s) for tests of the 
pQCD calculation and various lowest order cross sections, 
i.e. the measurement of i?had(s), in this energy region, 
the overlapping effects have to be clarified clearly. Near 
the DD threshold region, say above the ip^S) resonance 
and below the DD threshold, no direct measurement of 
the Ruds is currently available. The measurements of the 
R(s) [6(] in this region and in the region between 2.0 and 
3.55 GeV which were obtained from analyzing the 
data taken with the BES-II detector in the period from 
1998 to 1999, can not be used directly to test the validity, 
since those give the -Rhad(s) rather than the i? u ds- 

In this Letter, we report direct measurements of 
i? u ds near the DD threshold region, measurements of 
-^uds(c )+i/>(377o)( s ) above the DD threshold and the 
-Rhad(s) values at 68 energy points in the region from 
3.650 to 3.872 GeV. We here define the -R u ds(c)+i>(377o)( s ) 
to be the R(s) ratio accounting for the contributions from 
both the continuum hadron production and the decays 
for "0(3770) — > hadrons. Combining the measurements of 
the -R u ds(*j+v>(3770)( s ) an d the cross section for DD pro- 
duction [7j would give us some information about non- 
DD decays of V>(3770) [HE^ and improve our knowl- 
edge about the nature of V> (3770). The data samples used 
in the analysis were taken with the BES-II detector 11 1 
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FIG. 1: The distribution of the event vertex in z. 



The of!*, is determined by 
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where -/V^J is the number of the observed hadronic 
events, L is the integrated luminosity, Chad is the effi- 
ciency for the detection of inclusive hadronic events and 
e had ^ s t ne trigger efficiency for collecting hadronic events 
in online data acquisition system. 

The hadronic events are required to have more than 
2 good charged tracks satisfying selection criteria as de- 
scribed in Ref. [§] . To separate some beam-gas associated 
background events and cosmic rays background events 
from hadronic events we calculated the event vertex in 
the beam line direction. Fig. Q] shows the distribution of 
the event vertex of the accepted events. Using a Gaus- 
sian function to describe the hadronic events plus a sec- 
ond order polynomial for the background to fit the dis- 
tribution, we obtain the number, of the candidates 
for the hadronic events. This number of candidates for 
hadronic events contains contaminations from some back- 
ground sources such as e + e~ — > t + t~, e + e _ — > ( r y)e + e~ , 
e + e~ — > (7)/i + /i~, and two-photon processes. The num- 
ber of background events, AT b , due to these processes 
can be estimated by means of a Monte Carlo simula- 
tion Subtracting from yields the number 
of the observed hadronic events, The systematic 
uncertainty in measuring the produced hadronic events 
due to the hadronic event selection criteria is estimated 
to be - 2.5% i). 

The integrated luminosities of the data sets are de- 
termined using large-angle Bhabha scattering events as 
described in Ref. [9J]. The systematic uncertainty in the 
measured luminosities is estimated to be ~ 1.9% 0. 

The detection efficiency for hadronic events is deter- 
mined via a special Monte Carlo generator [l2T | in which 
the radiative corrections to a 2 order are taken into ac- 
count. These generated events are simulated with the 
GEANT3-based Monte Carlo simulation. The system- 
atic uncertainty in the efficiencies due to the generator is 
estimated to be ~ 2.0% for reconstruction of the inclu- 
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FIG. 2: (a) The efficiency versus the nominal cm. energy; 
(b) The ISR factor versus the nominal cm. energy (see text). 
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FIG. 3: The -R uds (c)+v<(3770) (s) or R heid (s) versus the cm. en- 
ergy (see text). 



and ip(2S) decays, and to be ~ 0.7% for reconstruction 
of the continuum hadronic events @. Fig. EJa) shows 
the Monte Carlo efficiencies for the detection of hadronic 
events produced at the different nominal cm. energies, 
where the error is statistical. 

The trigger efficiencies are measured to be e tl i g = 
(100.01-°;°)% for both the e+e" 



(7)e + e and e + e 



-0.5J 

hadrons events [9( . 

To get the crj| ad (s) the cgk|(s) has to be corrected for 
the initial state radiative and vacuum polarization cor- 
rections. There is no correction for final state radiation 
needed, since the crg ad (s) includes the final state photon 
radiation. The correction factor, (1 + 5(s)), is given by 



(1 + S(s)) 



q cxp (s) 
<T°( S ) ' 



where er cxp (s) is the expected cross section and cr°(s) 



production. The cr exp (s) can be written as 



- oxp (s) - 



dx ■ 



a°(s(l-x)) 
|l-H( S (l-x))P 



;F(x, S ), 



(4) 

where a°(s(l—x)) is the total lowest order cross section in 
the energy range from 0.28 GeV to y/s, F(x, s) is a sam- 
pling function [13| and | 1 _n( s ( 1 i- a ;))P ^ s ^ ne correction fac- 
tor for the effects of vacuum polarization including both 
the leptonic and hadronic terms in QED [tj [l3j]. As- 
suming that there are no other new structure and effects 
except the 0(3770) and continuum hadrons in the energy 
region from 3.729 to 3.872 GeV, we calculate the radia- 
tive correction factors. In the calculation of the cross 
sections, the 0(3770) and ip(2S) resonance parameters 
measured by the BES Collaboration jT5| are used. In- 



serting the resonance parameters of other 1 states [14 1 
and the i? uds = 2.218 ± 0.091 measured by the BES Col- 
laboration jj| near 3.65 GeV in Eqs. (6)-(25) of Ref. [§] 
as the initial input, we obtain the radiative correction 
factors at the 68 energy points as shown in Fig. &b)- 
The uncertainty in (1 + S(s)) is less than 1.5% [9j. 

The sum of the lowest order cross sections for contin- 
uum hadron and "0(3770) production is given by 



; uds(c)+V(3770) 



<d(s) 
(l + 6(s)) 



i(s), (5) 



in which XL^^Res i( s ) * s the contribution from the 1 
resonances [9j near the cm. energy of 3.70 GeV except 
0(3770). Eq. (5) indicates 



7 uds(c)+^(3770)( S ) — (^uds(c) +^V(3770)) 



■00, (6) 



where R 



V>(3770) 



(s) is the R(s) value due to the decays 



for 0(3770) -> hadrons. Dividing the v° ds{c)+4 , {3770) (s) 
b y a l+^-( s ) yields i?uds(c)+i/-(3770)(s) values as summa- 
rized in table [H where i? u ds(c)+v>(3770)( s ) means i? u ds(s) 
or ^uds(c)( s ) + ^(3770) ( s )j below the DD threshold 
-^uds( c)+i/>(377o)( s ) equals i? u d s (s); above the threshold 

#uds(c)-H/>(3770)(s) = ^uds(c)(s) + #V(3770) (s) IS the R(s) 

value including both the continuum hadron production 
and the decays of -0(3770) — > hadrons. As mentioned 
above, calculating the effects of the vacuum polarization 
on a(Af§) and a^ M needs the i?had(s) values, which can 
be obtained by dividing 



o-had(s) 



_obs 
'had 



(s) 



(1+5(8)) 



(7) 



by 



Table Q] also lists the i?had(s) values. 
(3) The first errors in the measured i? u ds(c)+i/>(3770) (s) and 
i?had(s) values listed in the table are statistical including 
the point-to-point systematic uncertainty arising from 
the statistical uncertainty in the measured luminosity 
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the second ones are common systematic, arising from 
the uncertainties in luminosity (~ 1.9%), in selection 
of hadronic event (~ 2.5%), in Monte Carlo Modeling 
(~ 2.0%), in radiative correction (~ 1.5%) for measure- 
ments of R(s) values off the 0(3770) resonance and in 
■0(3770) parameters (~ 2.7 %) for those in the 0(3770) 
resonance region between 3.74 to 3.82 GeV. Adding these 
uncertainties in quadrature yields the total systematic 
uncertainties to be ~ 4.0% and ~ 4.9% for these outside 
the 0(3770) resonance and within 0(3770) resonance re- 
gions, respectively. 

Weighting the first 8 -R u ds(c)+v>(3770)( s ) values below 
the DD threshold with the statistical error yields 

R uds = 2.141 ± 0.025 ±0.085, 

which can directly be compared to -Rud s CD = 2.15 ± 0.03 
calculated by pQCD @,[l|. 

Figure [3] displays the fl u d s (c)+^(3770)( s ) and/or R ha d(s) 
values from this measurements, together with those mea- 
sured by the MARK-I Collaboration [TJ], by the BES 
Collaboration from analysis of different data samples Q 
and the continuum -R^ 00 calculated by pQCD 0, , 
which is shown by two straight lines indicating ±lcr er- 
ror interval. There is one thing which we would like 
to point out. In the measurements of the R(s) val- 
ues, the measurements reported in Ref. Q are the quan- 
tity i?uds(c)+v(377o) (s) which equals i? uds below the DD 
threshold and equals i? u ds(c)(s) + Rip(377o)( s ) above the 
DD threshold; while the measurements reported from 
this work are the quantities i? u ds, -R u ds(c)+v>(377o)( s ) and 
-Rhad(s), the other measurements are only the quantity 

-Rhad(s). 

In summary, we measured the i? u ds(c)+v>(3770) ( s ) and 
-Rhad(s) values in the energy region from 3.650 to 3.872 
GeV with the systematic uncertainty of ~ 4.0% for the 
continuum i? u d s and the -Rhad(s) outside the 0(3770) res- 
onance regions, and with the uncertainty of ~ 4.9% for 
the -R u ds(c)+</.(3770)(s) and the i? ha d(s) in the 0(3770) res- 
onance region. These are improved measurements with 
respect to the earlier measurements @, [l6| in both the 
accuracy and the number of energy points. These im- 
proved measurements of the Rhad(s) are expected to pro- 
vide an improvement in the precision of the calculated 
values of a(M\) and a^ M [| [TJ, [Tj|. We obtained 
the continuum i? u d s near the DD threshold region to 
be R uds = 2.141 ± 0.025 ± 0.085, which is in excellent 
agreement with R^ D = 2.15 ± 0.03 0, HU predicted 
by pQCD in this energy region. This measured i? u d s can 
directly be used to evaluate the strong coupling constant 
a s (s) at the energy scale of ~ 3 GeV. In this analysis, we 
clarified the overlapping effects of different processes to 
the i?had(s). This is important for test of the validity of 
the pQCD calculation near the DD threshold region and 



for correctly using the data in calculation of a(Af§) and 
a SM = ( g - 2)/2 with the help of the pQCD prediction. 
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TABLE I: Summary of the i? u d s (c)+i/ , (3770) ( s ) an d -Rhad(s) values measured at 68 energy points. 



cm. energy ^ uds (c)+^(3770) (f) ^had (s) cm. energy fi uds (c)+^(3770) (s) #had Q) 
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± 


0, 


116 


2 


,380 


± 


0, 


106 


± 


0. 


.116 


3 


,7804 


3, 


,382 


± 


0. 


137 


± 


0. 


166 


3 


.385 


± 


0. 


137 


± 


0. 


166 


3 


.7500 


2 


,517 


± 


0. 


085 


± 


0, 


123 


2 


,525 


± 


0. 


.085 


± 


0. 


,123 


3 


,7808 


3 


336 


± 


0. 


129 


± 


0. 


163 


3 


.340 


± 


0. 


129 


± 


0. 


163 


3 


.7512 


2 


,637 


± 


0. 


,090 


± 


0, 


129 


2 


,644 


± 


0. 


,090 


± 


0. 


,129 


3 


,7810 


3 


,464 


± 


0. 


138 


± 


0. 


170 


3 


.468 


± 


0. 


138 


± 


0. 


,170 


3 


.7524 


2 


.615 


± 


0, 


.095 


± 


0, 


128 


2 


,622 


± 


0, 


.095 


± 


0, 


.128 


3 


.7812 


3 


396 


± 


0. 


130 


± 


0. 


166 


3 


.399 


± 


0. 


130 


± 


0, 


166 


3 


.7536 


2 


,652 


± 


0, 


.093 


± 


0, 


130 


2 


,659 


± 


0. 


.093 


± 


0. 


.130 


3 


.7814 


3 


,514 


± 


0. 


124 


± 


0. 


172 


3 


.518 


± 


0. 


124 


± 


0, 


.172 


3 


.7548 


2 


,733 


± 


0. 


093 


± 


0, 


134 


2 


,739 


± 


0. 


093 


± 


0, 


,134 


3 


.7816 


2 


,944 


± 


0. 


137 


± 


0. 


,144 


2 


.947 


± 


0. 


137 


± 


0. 


144 


3 


.7560 


2 


,585 


± 


0. 


090 


± 


0, 


127 


2 


.591 


± 


0. 


090 


± 


0. 


127 


3 


.7818 


3 


,140 


± 


0. 


125 


± 


0. 


,154 


3 


.143 


± 


0. 


125 


± 


0. 


154 


3 


.7572 


2 


,942 


± 


0. 


,107 


± 


0, 


144 


2 


.948 


± 


0. 


.107 


± 


0. 


,144 


3 


.7822 


3 


253 


± 


0. 


124 


± 


0. 


159 


3 


.257 


± 


0. 


124 


± 


0. 


,159 


3 


.7584 


3 


,025 


± 


0. 


108 


± 


0, 


148 


3 


.031 


± 


0, 


108 


± 


0. 


148 


3 


.7826 


3, 


,326 


± 


0. 


115 


± 


0. 


163 


3 


.329 


± 


0. 


115 


± 


0, 


163 


3 


.7596 


3 


,076 


± 


0. 


102 


± 


0, 


151 


3 


.082 


± 


0. 


102 


± 


0. 


151 


3 


.7838 


3 


,154 


± 


0. 


114 


± 


0. 


155 


3 


.157 


± 


0. 


114 


± 


0. 


155 


3 


.7608 


3 


,138 


± 


0. 


089 


± 


0, 


154 


3 


.143 


± 


0. 


,089 


± 


0. 


,154 


3 


.7850 


2 


,879 


± 


0, 


107 


± 


0. 


141 


2 


.882 


± 


0. 


107 


± 


0. 


,141 


3 


.7620 


2. 


,992 


± 


0. 


110 


± 


0, 


147 


2 


.998 


± 


0. 


,110 


± 


0. 


,147 


3 


.7862 


2 


,902 


± 


0, 


105 


± 


0. 


142 


2 


.905 


± 


0. 


105 


± 


0. 


,142 


3 


.7622 


3 


,207 


± 


0. 


114 


± 


0. 


157 


3 


.213 


± 


0. 


,114 


± 


0. 


,157 


3 


.7874 


2 


,957 


± 


0, 


111 


± 


0. 


145 


2 


.960 


± 


0. 


111 


± 


0. 


,145 


3 


.7634 


3 


,345 


± 


0. 


.122 


± 


0. 


164 


3 


.350 


± 


0. 


.122 


± 


0. 


.164 


3 


.7886 


2, 


,571 


± 


0, 


097 


± 


0. 


126 


2 


.574 


± 


0. 


097 


± 


0, 


.126 


3 


.7646 


3 


,585 


± 


0. 


.126 


± 


0. 


176 


3 


.590 


± 


0. 


.126 


± 


0. 


.176 


3 


.7898 


2, 


,576 


± 


0, 


099 


± 


0. 


126 


2 


.579 


± 


0. 


099 


± 


0, 


.126 


3 


.7658 


3 


.381 


± 


0. 


119 


± 


0. 


166 


3 


.386 


± 


0. 


,119 


± 


0. 


166 


3 


.7900 


2 


,849 


± 


0, 


106 


± 


0. 


,140 


2 


.852 


± 


0. 


106 


± 


0. 


140 


3 


.7670 


3 


,760 


± 


0. 


130 


± 


0. 


184 


3 


.764 


± 


0. 


130 


± 


0. 


,184 


3 


.7950 


2 


,751 


± 


0, 


101 


± 


0. 


135 


2 


.754 


± 


0. 


101 


± 


0. 


135 


3 


.7682 


3 


,451 


± 


0. 


124 


± 


0. 


169 


3 


.455 


± 


0. 


124 


± 


0. 


169 


3 


.8000 


2 


,212 


± 


0, 


091 


± 


0. 


108 


2 


.215 


± 


0. 


091 


± 


0. 


108 


3 


.7694 


3. 


,611 


± 


0. 


125 


± 


0. 


177 


3 


.615 


± 


0. 


,125 


± 


0. 


.177 


3 


.8100 


2 


,171 


± 


0, 


092 


± 


0. 


.087 


2 


.173 


± 


0. 


092 


± 


0, 


087 


3 


.7706 


3 


,580 


± 


0. 


123 


± 


0. 


175 


3 


.584 


± 


0. 


123 


± 


0. 


175 


3 


.8200 


2, 


367 


± 


0. 


109 


± 


0. 


.095 


2 


.369 


± 


0. 


109 


± 


0. 


095 


3 


.7714 


3 


,538 


± 


0. 


139 


± 


0. 


173 


3 


.543 


± 


0. 


139 


± 


0. 


,173 


3 


.8300 


2 


,354 


± 


0, 


101 


± 


0. 


094 


2 


.355 


± 


0. 


101 


± 


0. 


,094 


3 


.7716 


3 


,634 


± 


0. 


,146 


± 


0. 


178 


3 


.638 


± 


0. 


,146 


± 


0. 


,178 


3 


.8400 


2, 


,296 


± 


0, 


104 


± 


0. 


092 


2 


.297 


± 


0. 


104 


± 


0. 


,092 


3 


.7718 


3 


,939 


± 


0, 


.133 


± 


0, 


193 


3 


.943 


± 


0. 


.133 


± 


0. 


.193 


3 


.8500 


2, 


.372 


± 


0, 


115 


± 


0, 


095 


2 


.373 


± 


0. 


115 


± 


0, 


.095 


3 


.7720 


3 


636 


± 


0. 


.134 


± 


0, 


178 


3 


.640 


± 


0. 


.134 


± 


0. 


.178 


3 


.8600 


2 


.371 


± 


0, 


105 


± 


0, 


095 


2 


.372 


± 


0. 


105 


± 


0, 


095 


3 


.7722 


3 


,652 


± 


0, 


143 


± 


0, 


179 


3 


.656 


± 


0. 


143 


± 


0. 


.179 


3 


,8720 


2 


,308 


± 


0, 


117 


± 


0, 


092 


2 


.309 


± 


0. 


117 


± 


0. 


.092 



